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Abstract 

We present a computational study of the pump-probe response of a single atom to assess any 
microscopic nonlinear optical anisotropy due to freed electrons for non-resonant optical excitation. 
Using simulations of the Schrodinger equation for an atom exposed to a strong, short-duration 
linearly polarized pump, we calculate the induced dipole moment along probe directions parallel 
and perpendicular to the pump polarization. Our simulations show birefringence ratios of approx- 
imately 0.78-0.9 on the timescale of tens of femtoseconds following the excitation pulse. 

PACS numbers: 32.80.Rm,42.50.Hz 
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I. INTRODUCTION 



Isotropic dielectric media, such as gases or cubic crystals, can exhibit nonlinear birefrin- 
gence for incident field frequencies in the transparency region [1]. This nonlinear optical 
anisotropy (NOA) is revealed in pump-probe schemes in which a linearly polarized pump 
field is incident on the medium along with a co-propagating weak probe field polarized ei- 
ther parallel or perpendicular to the pump. For the case of a Kerr nonlinear medium, this 
yields a difference in refractive index experienced by the orthogonally polarized probes that 
is proportional to pump intensity. In other words, the orthogonally polarized probes induce 
different dipole moments along their respective directions in the presence of the pump. 

Various pump-probe schemes are at the heart of recent approaches to measure the non- 
linear opticalproperties and filamentation dynamics of atmospheric pressure gases, such as 
air or argon |2j-l5|. If the pump pulse is sufficiently intense to generate freed electrons via 
multi-photon or tunneling ionization, this raises the issue of whether these electrons also 
contribute to the NOA. Here, we use the term freed electrons to describe those that are 
no longer bound, but may still be correlated with the parent ion, as opposed to strictly 
free electrons. Physically, the linearly polarized pump pulse drives the freed electrons along 
the polarization direction, thereby creating an anisotropic momentum and spatial distribu- 
tion for the electrons. However, it is a common assumption that freed electrons exhibit 
an isotropic optical response. This is true if the freed electrons and associated ions have 
had ample time to undergo Coulomb scattering and produce a collective plasma whose mo- 
mentum distribution is isotropic. From simulations of the quantum Boltzmann equation 
for an atomic gas, it has been established js] that anisotropic features of the momentum 
distribution survive for a time of order t a ~ 1/oj p i, given by the inverse plasma frequency. 
A collective isotropic plasma appears for longer times t cp , with t a ~ 500 fs and t cp a few pi- 
coseconds for conditions typical of femtosecond filamentation in atmospheric pressure gases. 
Thus, for pulse durations greater than t cp , freed electrons should yield an isotropic optical 
response. However, for pulse durations less than t a , we expect that any anisotropy due to 
freed electrons will persist and may lead to a NOA. Photo-induced birefringence in semicon- 
ductors due to single-photon ionization provides an example of a NOA due to freed electrons 
for optical fields [7]. 

It has been recently shown [8| that an anisotropic electron momentum distribution can 
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reveal itself using a THz probe. In this paper, we present a computational study of the pump- 
probe response of a single atom, in order to assess any microscopic NOA due to freed electrons 
for optical fields. Essentially, we inquire if the freed electrons retain a "memory" of the linear 
polarization of the pump pulse that generated them. We restrict simulated pulse durations 
to less than t a , so that anisotropy of the electron distribution due to the pump pulse will be 
minimally degraded by electron-electron interactions, and so that a single atom approach is 
sufficient. In particular, using simulations of the single particle Schrodinger equation, the 
induced dipole moment along two probe directions is calculated to assess any NOA. Two 
pump models are examined: a delta-function pump of zero temporal duration, and a realistic 
pump with non-zero duration. In both cases, we find that a sizeable anisotropic nonlinear 
atomic response can occur due to freed electrons. Our hope is that these computational 
results provide impetus for future experimental and theoretical studies of NOA due to freed 
electrons on ultrashort time scales. In addition to being of fundamental interest, our results 
may have implications for the interpretation of ultrashort pump-probe experiments aimed 
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at measuring nonlinear optical properties in the presence of freed electrons (e.g. 
they rely on the assumption of an isotropic freed-electron response. Similarly, "plasma" 
isotropy is invariably assumed in computer aided modeling of femtosecond filamentation 
so the current model may require corrections. 

II. NUMERICAL MODEL 

Our starting point is the three-dimensional (3D) Schrodinger equation in atomic units 
(a.u.) for the wave function \P(r*, t) of an electron under the actions of the Coulomb potential 
due to the parent ion and the pump and probe fields 



-V-- + F(r,t) 
2 r 



% (1) 



where V(f,t) = f ■ [Ep(t) + E pr (t)] = \Vp + V pr \ is the interaction energy between the 
atom and pump (P) and probe (pr) fields in the dipole approximation. In the absence of 
pump and probe fields, the spherically symmetric ground state wave function is ty g (r,t) = 
ip g (r) exp(— ie g t), with if} g (r) = AAexp(— r), Af a normalization constant, and e g the ground 
state energy. 

We take the pump field to be linearly polarized, Ep{t) = e z E P (t), with e z the unit vector 
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along the z-axis. Then, employing cylindrical coordinates r = (p,8,z), the cylindrically 
symmetric solution ty(r,t) = ^p(p, z,t) of Eq. ([1]) in the presence of the pump beam alone 
evolves as id^ p/dt = Hp^> p with the Hamiltonian 

1 A 9 d d 2 \ 1 
H P = -n\ -ttP^t + 7TT - - + ( 2 ) 
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The probe is assumed much weaker than the pump, which makes \l/p the zeroth-order wave 
function, and the atom-probe interaction is treated to first-order in perturbation theory. The 
probe field is written as E pr {t) = e m E pr (t), where m = corresponds to the probe polarized 
along the z-axis, e*o = e z , and m = 1 to the probe polarized along the x-axis, e\ = e x . 
We express V pr = £, m cos(m9)E pr (t), with £ = z and = p. Writing the wave function 
as \I/(r, t) = typ(p,z,t) + ip m (p, z, t) cos(m6 l ), the Schrodinger equation for the first-order 
correction due to atom-probe interactions is 

.<9Vv, 

with the initial condition ip m (p, z, 0) = 0. We have built into this equation, the same probe 
field E pr (t) applied along both axes. 

Our approach involves generating numerical solutions for typ and ipo i, which are used to 
calculate the light-induced atomic dipole moments. In particular, the dipole moment due 
to the pump alone 

roc poo 

d P (t) = 2ir dpp / dz^fpz^fp, (4) 

JO J-oc 

with the resulting dipole moment directed along the z-axis. The probe-induced changes to 
the dipole moment 



ip m + £ m E pr (t)^ P (p, z,t), (3) 
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where d m {t) is directed along the z-axis (x-axis) for m = (m — 1) and corresponds to the 
raw data that we extract from the simulations. Furthermore, we also calculate the probe- 
induced changes to the dipole moment in the absence of the pump beam di^ (t) , which are 
obtained using Eq. (jHJ) with ^p — > The nonlinear change in the dipole moment is then 
given by 8d m (t) = [d m (t) — dm\t)}. We shall discuss further processing of the data as we 
proceed. 

Turning now briefly to some details of the numerical approach, the Schrodinger equations 
for ^ p and ipQ i are converted to systems of ordinary differential equations on a discrete grid 




100 200 300 400 500 600 700 
Time (a. u.) 



FIG. 1: (Color online) Scaled probe-induced dipole moment (dm (t) / Ep r>max ) versus time, for the 
cases (black crosses) m = (z-polarized probe) and (red line) m = 1 (x-polarized probe). This 
simulation verifies that isotropy is properly preserved in absence of the pump. 

in the two spatial coordinates (p, z). The Coulomb potential is discretized by asymptotic be- 
havior correspondence to ensure accuracy and improve computational efficiency [10!], and we 
begin the simulation in the ground state of the discretized system of equations. We perform 
a convergence study of measured observables to ensure that both radial and longitudinal grid 
sizes are sufficient to accommodate the expanding, pump-driven wave function. The results 
shown here use a grid resolution of 0.3 a.u. with a computational box size up to 9000 a.u 
by 250 a.u. in the z and p directions, respectively. Because we calculate and compare the 
induced dipole moments in different directions, our approach here is to avoid using absorbing 
boundaries which could affect results. Rather, we choose to restrict our simulations to times 
short enough such that a convergence study shows absence of boundary-related artifacts. 



III. RESULTS 



The goal of these calculations is to detect any anisotropy that may exist in the nonlinear 
response caused by a linearly polarized excitation. Thus, we verify that the model presented 
properly preserves isotropic properties where appropriate. To test the utility of the pertur- 
bative equations (J3J), we solve for the probe-induced dipole moments in the absence of the 
pump beam for off-resonant probe fields. Figure [T] shows the scaled, probe-induced dipole 
moments (dm (t) / E pr ^ max ) versus time for the case m = (z-polarized) and the case m = 1 
(x-polarized). A probe field E pr (t) is used with maximum amplitude E pr ^ max and a temporal 
profile matching d!"m (t) shown in the figure. We see that d ( 9 \t) = d[ 9 \t), meaning the linear 
probe response is isotropic, as expected. Furthermore, we find that dm\t) = aE pr (t) } with 
a > 4.5, very close to the analytic result for the off-resonant linear polarizability of the 
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FIG. 2: (Color online) Calculated probe-induced dipole moments (red line) do(t) and (black crosses) 
di(t) for a delta-function pump with (a) F = 0.041, and (b) F = 0.17. The probe has a center 
wavelength of X pr = 800 nm and has a temporal profile profile shown by the blue line in Fig. [3|a). 

atom. Thus, we have validation of the approach for the linear optical response. 

We now turn to simulations of the pump-probe response. The first model considered is for 
a delta-function pump pulse of the form Ep{t) = FS(t). Despite the extreme nature of this 
pulse, namely that it has zero temporal duration and hence contains all frequencies in equal 
proportion, we study this model for three main reasons. First, following the ^-polarized 
oump pulse, an atom initially in its ground state becomes \&p(r, t = 0+) = ipg{ r ) ex P( — iFz) 



ll| , which eliminates the need for a numerical solution to model the pump modification of 
the wave function. Second, this provides a simple model for the case when the pump and 
probe are temporally separated. In this limit, the probe experiences the atomic state left 
behind by the pump pulse but with no NOA due to the Kerr effect, since the pump has 
passed. Finally, we are optimistic that the simplifications this model affords will eventually 
yield analytic insights into the physics underlying the NOA. 

Figure [2] shows the calculated raw data in the form of the probe-induced dipole moments 
do t i(t) for a delta-function pump pulse with F = 0.041 and F = 0.17. In each case, the probe 
pulse has a center wavelength of X„ r = 800 nm and a temporal profile shown in Fig. E^a). 



Based on the results from Ref. the ionization probability is 0.1% for F = 0.041 and 
1.8% for F = 0.17. In each case the pump leaves behind the probability of a freed electron 
that can modify the probe. We note that for the relatively weak pump of F — 0.041 in Fig. 
12(a), the temporal profile of the dipole moments c?o,i(t) closely follow the probe profile in 
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FIG. 3: (Color online) Filtered nonlinear changes in dipole moment for (red line) 5do(t) and (black 
crosses) 5d\(t) for a delta-function pump with (a) F = 0.041 and (b) F = 0.17. The blue dashed 
line in (a) shows the temporal profile of the probe in arb. units. 

Fig. El^a), indicating a weak nonlinear response, whereas the temporal profiles of the dipole 
moments do,x(t) for F = 0.17 in Fig. EJ^b) are clearly different and show a much stronger 
nonlinear effect. 

We are primarily interested in the nonlinear change in dipole moment for frequencies 
in the vicinity of the probe frequency u pr , and for this we employ the filtering procedure 



described in Ref. 12]. Obtaining the nonlinear change in dipole moment 5d m (u) from 5d m (t) 
via a Fourier transform, a band-pass filter J-(u) — u pr ) centered on u pr is applied to isolate 
the nonlinear response in the vicinity of the probe. The inverse Fourier transform yields the 
filtered nonlinear change in dipole moment Sd m (t). The robustness of this filtering procedure 
has been tested by considering various Gaussian and super-Gaussian filters. Figure [3] shows 
5do t i(t) for F = 0.041 and F = 0.17. There are several key features of these results. First, 
despite the differences in the raw data in Figs. EJ^a) and|2]^b), the filtered nonlinear changes 
in dipole moment in Fig. [3] are very similar in temporal profile, though they have different 
magnitudes. This highlights the importance of subtracting the probe-induced change in the 
absence of the pump. Second, 5d 0; i(t) are 180° out of phase with the probe pulse [see Fig. 
Efa)], which is consistent with an optical response due to freed electrons. Finally, for both 
examples in Fig. [HI Sdo(t) along the z-axis is seen to be greater than Sdi(t) along the x-axis, 
thereby demonstrating a NOA due to the freed electrons, a key result of this paper. There 
is no NOA due to the Kerr effect as the pump and probe do not overlap. 
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FIG. 4: Birefringence ratio B r = 5di :Pea k/Sdo tPea k as a function of the pump parameter F, with 
fidm^peak the peak value of 5d m {oj) near u = uj pr . 

To further explore the NOA revealed above, Fig. H] shows the birefringence ratio B r = 
Sdi iPea k/6do tPea i~ as a function of the pump parameter F, with 5d m ^ pea k the peak value of 
Sd m (u) near uj = u pr . We checked that the results in Fig. @] are largely unchanged by 
variations in the filtering procedure. The calculation is limited at small values of F (< 0.001) 
by numerical noise arising from the subtraction of the relatively large probe-induced dipole 
moment in the absence of the pump. At large values of F (> 0.5), the momentum supplied to 
the atomic wave function due to the pump pushes it strongly in one direction. Therefore, the 
large zero-frequency component dominates the spectrum, making isolation of the nonlinear 
change to the dipole moment in the vicinity of the probe frequency increasingly difficult. 
Over the range shown of F, the birefringence ratio varies between B r « 0.78 — 0.84. To 
put this in context, using the results from Ref. Ill ], the ionization probability is 0.1% 
for F = 0.041 and 1.8% for F = 0.17, so that these significant birefringence ratios arise for 
ionization probabilities that are not particularly large and are compatible with the conditions 
present during laser filamentation. 

We next consider a more realistic model with a pump pulse of center wavelength Xp = 800 
nm, peak intensity 10 14 W/cm 2 , and 13 fs duration along with a weak probe pulse of center 
wavelength X pr = 400 nm and duration 30 fs. We choose the probe at the second-harmonic 
of the pump to aid in filtering out the change in dipole moment in the vicinity of the probe 
wavelength. The envelopes of the pump field and scaled probe field are shown in Fig. 0(a), 
which are seen to overlap. We also note that the pulse durations are much shorter than the 
time scale t a ~ 500 fs, on which anisotropic features in the electron momentum distribution 
due to the pump pulse are expected to disappear. 

Figure E](a) shows the filtered nonlinear changes in dipole moment Sd m (t) versus time 
with the probe field directed along the z-axis (a;- axis) for to = (to = 1). For t > 1200 a.u., 
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FIG. 5: (Color online) Filtered nonlinear changes in dipole moment for (red line) Sdo(t) and (black 
crosses) 8di(t). Also shown are the envelopes of the (blue dashed line) scaled-probe field and (green 
dotted line) pump field with (a) intensity 10 14 W/cm 2 and duration 13 fs, and (b) intensity 2- 10 14 
W/cm 2 and duration 8 fs. 

after the pump pulse has passed, we again see that 5d (t) > 5di(t), indicating the presence 
of a NOA. Furthermore, the filtered nonlinear changes in dipole moment at these times are 
180° out of phase with the probe field so that the NOA is due to freed electrons. Inspection 
of the results for t > 1200 a.u. yields a birefringence ratio B r ~ 0.88, closer to unity than 
for the delta-function excitation, but still significant. For t < 750 a.u., it appears that the 
NOA is much more prominent with B r « 1/3. During this time region the filtered nonlinear 
changes in dipole moment are in-phase with the probe field and the probability of ionization 
is small, indicating the observed NOA is dominantly due to the nonlinear Kerr effect 13]. 
In the intermediate time region 750 < t < 1200 a.u., the NOA is due to a combination of 
the Kerr effect and freed electrons. 

In Fig. El(b), the pump pulse duration is shortened to 8 fs and the intensity increased 
to 2 • 10 14 W/cm 2 . The higher intensity induces a NOA from nearly the beginning of 
the simulation that is a combination of the Kerr effect and freed electrons, so that the 
birefringence due to the Kerr effect alone is no longer dominant. After the pump pulse has 
passed, a NOA due to freed electrons alone is observed again, but for times t > 1500 a.u. 
the response becomes isotropic. Averaging over 1200 < t < 2300 a.u., B r ~ 0.92. This 
behavior is consistent with the results of Fig. HJ where the birefringence weakens for higher 
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pump values. 



IV. CONCLUSION 

In summary, we have presented computational results based on 3D time-domain 
Schrodinger equation simulations to demonstrate that the pump-probe response of a sin- 
gle atom displays a sizeable nonlinear optical anisotropy (NOA) due to freed electrons for 
parameters compatible with those in light filamentation in atmospheric gases. The origin of 
the effect is in the anisotropic momentum distribution and spatial correlation in electrons 
freed by a strong linearly polarized pump. We found a NOA clearly present for two pump 
models, namely a delta-function pulse, as well as a more realistic optical pump pulse of finite 
duration. Due to numerical limitations, we can only detect a NOA up to times of a few 
tens of femtosecond after excitation. However, it is conceivable that it can persist longer, 
being bounded from above by the inverse plasma frequency for a gas of atoms |f|. The 
NOA revealed here could play a major role in the dynamics of vector light filamentation as 
the freed electron contribution will modify the polarization of the propagating filament in a 
manner not previously considered. 
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